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Abstract — In this paper, we present robust joint non-linear 
transceiver designs for multiuser multiple-input multiple-output 
(MIMO) downlink in the presence of imperfections in the channel 
state information at the transmitter (CSIT). The base station 
(BS) is equipped with multiple transmit antennas, and each user 
terminal is equipped with one or more receive antennas. The 
BS employs Tomlinson-Harashima precoding (THP) for inter- 
user interference pre-cancellation at the transmitter. We consider 
robust transceiver designs that jointly optimize the transmit 
THP filters and receive filter for two models of CSIT errors. 
The first model is a stochastic error (SE) model, where the 
CSIT error is Gaussian-distributed. This model is applicable 
when the CSIT error is dominated by channel estimation error. 
In this case, the proposed robust transceiver design seeks to 
minimize a stochastic function of the sum mean square error 
(SMSE) under a constraint on the total BS transmit power. 
We propose an iterative algorithm to solve this problem. The 
other model we consider is a norm-bounded error (NBE) model, 
where the CSIT error can be specified by an uncertainty set. 
This model is applicable when the CSIT error is dominated 
by quantization errors. In this case, we consider a worst-case 
design. For this model, we consider robust i) minimum SMSE, ii) 
MSE-constrained, and Hi) MSE-balancing transceiver designs. 
We propose iterative algorithms to solve these problems, wherein 
each iteration involves a pair of semi-definite programs (SDP). 
Further, we consider an extension of the proposed algorithm 
to the case with per-antenna power constraints. We evaluate 
the robustness of the proposed algorithms to imperfections in 
CSIT through simulation, and show that the proposed robust 
designs outperform non-robust designs as well as robust linear 
transceiver designs reported in the recent literature. 

Index Terms — Multiuser MIMO downlink, non-linear pre- 
coding, imperfect CSIT, Tomlinson-Harashima precoder, joint 
transceiver design. 

I. Introduction 

Multiuser multiple-input multiple-output (MIMO) wireless 
communication systems have attracted considerable interest 
due to their potential to offer the benefits of spatial diversity 
and increased capacity [1], [2]. Multiuser interference limits 
the performance of such multiuser systems. To realize the 
potential of such systems in practice, it is important to devise 
methods to reduce the multiuser interference. Transmit-side 
processing at the base station (BS) in the form of precoding 
has been studied widely as a means to reduce the multiuser 
interference [2]. Several studies on linear precoding and non- 
linear precoding (e.g., Tomlinson-Harashima precoder (THP)) 
have been reported in the literature [3], [4]. Joint design of 
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both transmit precoder and receive filter can result in im- 
proved performance. Transceiver designs that jointly optimize 
precoder/receive filters for multiuser MIMO downlink with 
different performance criteria have been widely reported in 
the literature [5]-[ll]. An important criterion that has been 
frequently used in such designs is the sum mean square error 
(SMSE) [6]-[9]. Iterative algorithms that minimize SMSE 
with a constraint on total BS transmit power are reported 
in [6], [7]. These algorithms are not guaranteed to converge 
to the global minimum. Minimum SMSE transceiver designs 
based on uplink-downlink duality, which are guaranteed to 
converge to the global minimum, have been proposed in [8], 
[9]. Non-linear transceivers, though more complex, result in 
improved performance compared to linear transceivers. Studies 
on non-linear THP transceiver design have been reported in 
the literature. An iterative THP transceiver design minimizing 
weighted SMSE has been reported in [10]. References [8] and 
[11], which primarily consider linear transceivers, present THP 
transceiver optimizations also as extensions. In [12], a THP 
transceiver design minimizing total BS transmit power under 
SINR constraints is reported. 

All the studies on transceiver designs mentioned above 
assume the availability of perfect channel state information 
at the transmitter (CSIT). However, in practice, the CSIT is 
usually imperfect due to different factors like estimation error, 
feedback delay, quantization, etc. The performance of precod- 
ing schemes is sensitive to such inaccuracies [13]. Hence, it 
is of interest to develop transceiver designs that are robust 
to errors in CSIT. Linear and non-linear transceiver designs 
that are robust to imperfect CSIT in multiuser multi-input 
single-output (MISO) downlink, where each user is equipped 
with only a single receive antenna, have been studied [14]- 
[19]. Recently, robust linear transceiver designs for multiuser 
MIMO downlink (i.e., each user is equipped with more than 
one receive antenna) based on the minimization of the total 
BS transmit power under individual user MSE constraints and 
MSE-balancing have been reported in [20]. However, robust 
transceiver designs for non-linear THP in multiuser MIMO 
with imperfect CSIT, to our knowledge, have not been reported 
so far, and this forms the main focus of this paper 

In this paper, we consider robust THP transceiver designs 
for multiuser MIMO downlink in the presence of imperfect 
CSIT. We consider two widely used models for the CSIT error 
[21], and propose robust THP transceiver designs suitable for 
these models. First, we consider a stochastic error (SE) model 
for the CSIT error, which is applicable in TDD systems where 
the error is mainly due to inaccurate channel estimation (in 
TDD, the channel gains on uplink and downlink are highly 
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Fig. 1. Multiuser MIMO downlink system model with Tomlinson-Harashima preceding. 



correlated, and so the estimated channel gains at the transmitter 
can be used for precoding purposes). The error in this model 
is assumed to follow a Gaussian distribution. In this case, 
we adopt a statistical approach, where the robust transceiver 
design is based on minimizing the SMSE averaged over the 
CSIT error To solve this problem, we propose an iterative 
algorithm, where each iteration involves solution of two sub- 
problems, one of which can be solved analytically and the 
other is formulated as a second order cone program (SOCP) 
that can be solved efficiently. Next, we consider a norm- 
bounded error (NBE) model for the CSIT error, where the error 
is specified in terms of uncertainty set of known size. This 
model is suitable for FDD systems where the errors are mainly 
due to quantization of the channel feedback information [17]. 
In this case, we adopt a min-max approach to the robust 
design, and propose an iterative algorithm which involves the 
solution of semi-definite programs (SDP). For the NBE model, 
we consider three design problems: i) robust minimum SMSE 
transceiver design ii) robust MSE-constrained transceiver de- 
sign, and in) robust MSE-balancing transceiver design. We 
also consider the extension of the robust designs to incorporate 
per-antenna power constraints. Simulation results show that the 
proposed algorithms are robust to imperfections in CSIT, and 
they perform better than non-robust designs as well as robust 
linear designs reported recently in the literature. 

The rest of the paper is organized as follows. The system 
model and the CSIT error models are presented in Section |ll] 
The proposed robust THP transceiver design for SE model 
of CSIT error is presented in Section Hill The proposed 
robust transceiver designs for NBE model of CSIT error are 
presented in Section |IV] Simulation results and performance 
comparisons are presented in Section [V] Conclusions are 
presented in Section [VT] 

II. System Model 

We consider a multiuser MIMO downlink, where a BS 
communicates with M users on the downlink. The BS employs 
Tomlinson-Harashima precoding for inter-user interference 
pre-cancellation (see the system model in Fig. [TJ. The BS 
employs Nt transmit antennas and the fcth user is equipped 



with N,.^ receive antennas, 1 < fc < M. Let Ufc denot^H 
the Lfc X 1 data symbol vector for the fcth user, where Lk, 
k = 1, 2, • • • , M, is the number of data streams for the fcth 
user Stacking the data vectors for all the users, we get the 
global data vector u — [uj , 



uIj] . The output of the fcth 



user's modulo operator at the transmitter is denoted by v^. 
Let Bfc € C^'^^'' represent the precoding matrix for the fcth 
user The global precoding matrix B = [Bi,B2,-- - ,Bj\/]. 
The transmit vector is given by 



= Bv, 



(1) 



where v = [vj, ■ ■ ■ , vj^]^. The feedback filters are given by 



Q.k,i ■ ■ 'iifc.fc-i 







, 1 < fc < M, (2) 



where Gj.^ 



••LkXLi 



, perform the interference pre- 
subtraction. We consider only inter-user interference pre- 
subtraction. When THP is used, both the transmitter and the 
receivers employ the modulo operator, Mod( ). For a com- 
plex number x, the modulo operator performs the following 
operation 



Mod(a;) 
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-ja 




1 


a 


^2 


a 


^2 



(3) 



where j = ^/—l, and a depends on the constellation [22]. For 
a vector argument x ~ [xi X2 ■ ■ ■ xn]'^, 



Mod(x) = [Mod(a;i) Mod(a;2) 
The vectors and are related as 



Mod(a;Ar)]^ . (4) 



Mod 



fe-i 



■k,3 ^3 



(5) 



The fcth component of the transmit vector x is transmitted 
from the fcth transmit antenna. Let denote the Nr^ x Nt 



We use the following notation: Vectors are denoted by boldface lowercase 
letters, and matrices are denoted by boldface uppercase letters. [.]^, [.]^, and 
denote transpose, Hermitian, and pseudo-inverse operations, respectively. 
[A.]ij denotes the element on the ith row and jth column of the matiix A. 
vec(.) operator stacks the columns of the input matrix into one column- vector. 
II ■ \\p denotes the Frobenius norm, and E{ } denotes expectation operator. 
A >: B implies A — B is positive semi-definite. 
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channel matrix of the fcth user. The overall channel matrix is 
given by 



H 



The received signal vectors are given by 

Yk = HfeBv + rife, l<k<M. 
The fcth user estimates its data vector as 



(6) 
(7) 



{CkYk) mod a 
= (CfeHfeBv + Cfcrifc) mod a, l<k<M, (8) 

where Cfc is the Lk x Nr^ dimensional receive filter of the fcth 
user, and is the zero-mean noise vector with E{nfen^} = 
CT^I. Stacking the estimated vectors of all users, the global 
estimate vector can be written as 

u = (CHBv + Cn) mod a, (9) 

where C is a block diagonal matrix with C/c, 1 < fc < Af 
on the diagonal, and n = [n^, • • • , n^^]^. The global receive 
matrix C has block diagonal structure as the receivers are 
non-cooperative. Neglecting the modulo loss, and assuming 
Ejvj.vl'^} — I, we can write MSE between the symbol vector 
Ufe and the estimate Ufe at the fcth user as [10] 

efc = E{||ufc-Ufe||2} 

(CfeHfeB — Gfc)(CfcHfeB — Gfe)^ + cr^CfeCf^ , 

1 <k<M, (10) 



where Gt 



G 



G 



■k,k-l 







A. CSIT Error Models 

We consider two models for the CSIT error. In both the 
models, the true channel matrix of the fcth user, H^, is 
represented as 

Hfc = Hfc + Efe, l<fc<M, (11) 

where is the CSIT of the fcth user, and is the CSIT 
error matrix. The overall channel matrix can be written as 



H 



H + E, 



(12) 



where H 



[Hf Hi 



and E 



[E^ E2 ■ ■ ■ E^^]^. In a stochastic error (SE) model, 
Efc is the channel estimation error matrix. The error matrix 
Efe is assumed to be Gaussian distributed with zero mean and 
EjEfcE^} — a'^lN^^N^i^- This statistical model is suitable 
for systems with uplink-downlink reciprocity. We use this 
model in Sec. |lll] An alternate error model is a norm-bounded 
error (NBE) model, where 



\Ek\\F < 4, l<fc<M, 



(13) 



or, equivalently, the true channel belongs to the uncertainty 
set TZk given by 

Uk - {C|C = Hfe + Efc, \\Ek\\F < Sk}, 1 < fc < A/, (14) 

where Sk is the CSIT uncertainty size. This model is suitable 
for systems where quantization of CSIT is involved [17]. We 



use this model in Sec. |IV[ 

III. Robust Transceiver Design with Stochastic 
CSIT Error 

In this section, we propose a transceiver design that min- 
imizes SMSE under a constraint on total BS transmit power 
and is robust in the presence of CSIT error, which is assumed 
to follow the SE model. This involves the joint design of the 
precoder B, feedback filter G, and receive filter C. When E, 
the CSIT error matrix, is a random matrix, the SMSE is a 
random variable. In such cases, where the objective function 
to be minimized is a random variable, we can consider the 
minimization of the expectation of the objective function. In 
the present problem, we adopt this approach. Further, the 
computation of the expectation of SMSE with respect to E 
is simplified as E follows Gaussian distribution. Following 
this approach, the robust transceiver design problem can be 
written as 



EE{smse} 
subject to Tr (BB^) < P„ 



mm 

B,C,G 



(15) 



where Pmax is the limit on the total BS transmit power, 
and minimization over B, C, G implies minimization over 
Bi,Ci,Gi. 1 < i < M. Incorporating the imperfect CSIT, 
H = H + E, in ([TOl), the SMSE can be written as 

smse = E{||u-u|p} 

tr[(Cfc(Hfc + Ek)B - Gfe)(Cfe(Hfe + Efe)B - G^ 



M 



k=l 



(16) 



Averaging the smse over E, we write the new objective 
function as 



A 



Eelsmse} 



M 

^tr[(CfcHfcB - Gk){CkUkB - Gk)" 



k=l 



(a|tr(BB«) + a2)CfcCf 



(17) 



Using the objective function fj,, the robust transceiver design 
problem can be written as 



mm 

B,C,G 



(18) 



Bjlj' < Pmax- 



subject to 

From ( [TtI i. we observe that fi is not jointly convex in B, G, 
and C. However, it is convex in B and G for a fixed value 
of C, and vice versa. So, we propose an iterative algorithm 
in order to solve the problem in ( fTSl ). where each iteration 
involves the solution of a sub-problem which either has an 
analytic solution or can be formulated as a convex optimization 
program. 

A. Robust Design of G and C Filters 

Here, we consider the design of robust feedback and receive 
filters, G and C, that minimizes the smse averaged over E. 
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For a given B and C^, as we can see from ( fTTj i. the optimum 
feedback filter Gf.j, 1 < k < M, j < k, is given by 

G,.,^. = CfeHfeB,. (19) 

Substituting the optimal Gj. ^ given above in ( [TtI i. the objec- 
tive function can be written as 



M 

E4( 



CfcHfcBfc — I)(CfcHfcBA 



fc=l 



+ ^ (CfcHfcBj)(CfeHfcBj)"^ 
+ (<7|tr(BB^)+a^)C,Cf . 



(20) 



In order to compute the optimum receive filter, we differentiate 
(|20] | with respect to Cfe, 1 < fc < M, and set the result to 
zero. We get 



M 



C,(h,( B,Bf)Hf 

j=k+i 



+ {al + al\\B\\' 
From the above equation, we get 



1 < fc < M. (21) 



A/ 

H.( Y: B,B-)H, 



<yUB\\l)I 



1 < fc < M. (22) 



We observe that the expression for the robust receive filter in 
(l22T l is similar to the standard MMSE receive filter, but with 
an additional factor that account for the CSIT error. In case 
of perfect CSIT, aE — and the expression in ( l22l i reduces 
to the MMSE receive filters in [10], [12]. 

B. Robust Design of B Filter 

Having designed the feedback and receive filter matrices, 
G and C, for a given precoder matrix B, we now present the 
design of the robust precoder matrix for given feedback and 
receive filter matrices. Towards this end, we express the robust 
transceiver design problem in ( fTSl l as 



mm 

b,c,g 



J2 W^khk - gkW + (4l|b|P + al)\\c,r (23) 



k=l 



subject to ||b||^ < P,nax, 

where = (B^ C^,), hj. = vec(Hfc), b = vec(B), Cj. = 
vec(Cfc), gk — vec(Gfe), and = vec(Hfe). Minimization 
over b, c, g denotes minimization over b^, c^, g^, 1 < i < M. 
For given C and G, the problem given above is a convex 
optimization problem. The robust precoder design problem, 
given C and G, car|^^e written as 

• f]||D,h,-g,|p + 4||b|p||cfe||2 

k=l 

(24) 



mm 

b 



-o--«l|cfc 



subject to 



< Prr 



As the last term in ( l24l i does not affect the optimum value of 
b, we drop this term. Dropping this term and introducing the 
dummy variables t^, J'fc, 1 < fc < M, the problem in ( |24] | can 
be formulated as the following convex optimization problem: 



mm 

subject to 



M 

E 



tk + cr£;||Cfc|| rfc 



(25) 



|Dfchfc-gfe||2 <tfc, 

l|bf <r,, 

Tk < Pmax, 1 < fc < M. 

The constraints in the above optimization problem are rotated 
second order cone constraints [23]. Convex optimization prob- 
lems like that in (|25] | can be efficiently solved using interior- 
point methods [23], [24]. 

C. Iterative Algorithm to Solve U5i 

Here, we present the proposed iterative algorithm for the 
minimization of the SMSE averaged over E under total BS 
transmit power constraint. In each iteration, the computations 
presented in subsections IIII-AI and IIII-BI are performed. In the 
(n + l)th iteration, the value of B, denoted by B"+^, is the 
solution to the following problem 

B"+i = argmin ^(B,C",G"), (26) 

B:Tr(BB'^)<P,„„^ 

which is solved in the previous subsection. Having computed 
■gn+i^ Qn+i jjjg sQjution to the following problem: 



= argmin/i(B"+\C,G"), (27) 
c 



and its solution is given in (l22l l. Having computed B"+^ and 
C"+^, G"+^ is the solution to the following problem: 



G 



ri+l _ 



= argmin/i(B"+\C"+\G) 

G 



(28) 



and its solution is given in (fT9] l. As the objective function 
in ([TtI i is monotonically decreasing after each iteration and 
is lower bounded, convergence is guaranteed. The iteration 
is terminated when the norm of the difference in the results 
of consecutive iterations are below a threshold or when the 
maximum number of iterations is reached. We note that the 
proposed algorithm is not guaranteed to converge to the global 
minimum. 

IV. Robust Transceiver Designs with 
Norm-Bounded CSIT Error 

When the receivers quantize the channel estimate and send 
the CSI to the transmitter through a low-rate feedback channel, 
we can model the error in CSI at the transmitter by the NBE 
model [17]. In such cases, it is appropriate to consider the 
min-max design, where the worst-case value of the objective 
function is minimized. In this section, we address robust 
transceiver designs in the presence of a norm-bounded CSIT 
error Specifically, we consider i) a robust SMSE transceiver 
design, ii) a robust MSE-constrained transceiver design, and 
iii) a robust MSE-balancing (min-max fairness) design. 
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A. Robust SMSE Transceiver Design 



Here, we consider a min-max design, wherein the design 
seeks to minimize the worst case SMSE under a total BS 
transmit power constraint. This problem can be written as 



mm 

B,C,G 



max smse(B,C,G,E) (29) 

||Efc||<5fc,Vfc 



subject to tr(BB'^) < Pmax- 

The above problem deals with the case where the true channel, 
unknown to the transmitter, may lie anywhere in the uncer- 
tainty region. In order to ensure, a priori, that MSE constraints 
are met for the actual channel, the precoder should be so 
designed that the constraints are met for all members of the 
uncertainty set. This, in effect, is a semi-infinite optimization 
problem [25], which in general is intractable. We show, in 
the following, that an appropriate transformation makes the 
problem in ( |29] l tractable. We note that the problem in 
can be written as 



mm 

b,c,g,t 



M 
k=l 



subject to ||Dfe(h;; 



(30) 

-ek)-gk\\^ +<jl\\ck\\^ <tk, (31) 
V||efc||<4, l<fc<M, 



llb||^ < P„,ax, 

where efc = vec(Efc). The first constraint in ( [30] l is convex 
in B and for a fixed value of and vice versa, but 
not jointly convex in B, Gfe and C^. Hence, to design the 
transceiver, we propose an iterative algorithm, wherein the 
optimization is performed alternately over {B, G} and {C}. 

7 ) Robust Design of B and G Filters: For the design of 
the precoder matrix B and the feedback filter G for a fixed 
value of C, the second term in the left hand side of the first 
constraint in ( [30] l is not relevant, and hence we drop this term. 
Invoking the Schur Complement Lemma [26], and dropping 
the second term, we can write the constraint in (|30] | as the 
following linear matrix inequality (LMI): 

H 



Dfc(hfe + efc) - gfe 



Dfc(hfc 



- efc) 
I 



gfc 



h 0. (32) 



Hence, the robust precoder and feedback filter design problem, 
for a given C, can be written as 



mm 

B.G,t 



subject to 



M 
k=l 

ife [Dfehfe - gfc] 

[Dfchfc - gfc] I 

V||efe||<4, l<k<M, 



(33) 



H 



h 0, 



||b|| < y/Prrtax, 

where hfc = hfe + Gfe. From ( [32] i, the first constraint in 
can be written as 



A >z P^XQ + Q^X^P, 



(34) 



where 



A = 



tk 

Dfehfe - 



Dfchfc - gfc 
I 



H 



(35) 



P = [0 Df], X = Gfe, and Q = - [1 0]. Having 
reformulated the constraint as in ( |34] i, we can invoke the 
following Lemma [27] to solve the problem in ( [33] l. 
Lemma 1: Given matrices P, Q, A with A = A^, 



A ^ P^XQ + Q^X^P, 
if and only if 3 A > such that 

A - AQ^Q -pV" 
-pP AI 



VX:||X||<p, (36) 



> 0. 



(37) 



□ 



Applying Lemma 1, we can formulate the robust precoder 
design problem as the following convex optimization problem: 



where 



Mfe 



mm 

B,G,t,,3 

subject to 



tk — Pk 

(Dfehfe - gfc) 




M 

k=l 

Mfe h 13k > 0, Vfc, 

||b|| < VPmax, 



(38) 



(Dfchfe - gfe 



-4Df 





-(^feDfe 
Pkl 



(39) 



2) Robust Design of Filter Matrix C: In the previous 
subsection, we considered the design of the B and G matrices 
for a fixed C. Here, we consider the robust design C for given 
B and G. This design problem can be written as 



mm 



M 

E 

fc=i 



tk 



(40) 



subject to 



||Dfe(hfe+Gfe)-gfe|| 
V||Efe||<4, l<k<M. 



crlWckW^ < tk, 



Applying the Schur Complement Lemma, we can represent 
the first constraint in ( l40l l as 



tk 

Dfe(hfe + Gfc) 

fnCfe 



Dfc(hfc + Gfe) 

fnCfc 



n H' 



gfe 



gfe 



h 0. (41) 



The second inequality in the above problem, like in the 
precoder design problem, represents an infinite number of 
constraints. To make the problem in ( |40] i tractable, we again 
invoke Lemma 1. Following the same procedure as in the 
precoder design, starting with dTIT i. we can reformulate the 
robust receive filter design as the following convex optimiza- 
tion problem: 



mm 
c.t.x 



subject to 



M 

Y^k 

fe=l 

Nfe h 0, 



(42) 



V/c, 
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where 







(Dfehfe - gfc 



H 



I 







,(43) 



where Fi. = 







5j Iterative Algorithm to Solve ( \29t : In the previous subsec- 
tions, we described the design of B and G for a given C, and 
vice versa. Here, we present the proposed iterative algorithm 
for the minimization of the SMSE under a constraint on the 
total BS transmit power, when the CSIT error follows NBE 
model. The algorithm alternates between the optimizations of 
the precoder/feedback filter and receive filter described in the 
previous subsections. At the (n + l)th iteration, the value 
of B, denoted by B"+^, is the solution to problem ( [38] l, 
and hence satisfies the BS transmit power constraint. Having 
computed B"+^, C"+^ is the solution to the problem in ( |42] |. 
So J(B"+i,C"+i) < J(B"+\C") < J(B",C"), where 



J(B,C) 



max {smse(B,C,G,E)}. (44) 



The monotonically decreasing nature of J(B",C"), together 
with the fact that J(B",C") is lower-bounded, implies that 
the proposed algorithm converges to a limit as n ^ oo. The 
iteration is terminated when the norm of the difference in the 
results of consecutive iterations are below a threshold or when 
the maximum number of iterations is reached. This algorithm 
is not guaranteed to converge to the global minimum. 

4) Transceiver Design with Per-Antenna Power Con- 
straints: As each antenna at the BS usually has its own am- 
plifier, it is important to consider transceiver design with con- 
straints on power transmitted from each antenna. A precoder 
design for multiuser MISO downlink with per-antenna power 
constraint with perfect CSIT was considered in [28]. Here, 
we incorporate per-antenna power constraint in the proposed 
robust transceiver design. For this, only the precoder matrix 
design dSSl l has to be modified by including the constraints on 
power transmitted from each antenna as given below; 



mm 

b 



subject to 



M 

fc=i 



(45) 



i(/.,,Br<Pfe, i<k<M, 



where 0^, = [Oixfe-i 
be computed using ( |42] |. 



1 Oi 



xNt 



The receive filter can 



B. Robust MSE-Constrained Transceiver Design 

Transceiver designs that satisfy QoS constraints are of inter- 
est. Such designs in the context of multiuser MISO downlink 
with perfect CSI have been reported in the literature [29], [30], 
[31]. Robust linear precoder designs for MISO downlink with 
SINR constraints are described in [32]. Here, we address the 



problem of robust THP transceiver design for multiuser MIMO 
with MSE constraints in the presence of CSI imperfections. 
THP designs are of interest because of their better performance 
compared to the linear designs. 

When the CSIT is perfect, the transceiver design under MSE 
constraints can be written as 



mm 

B,G,C 

subject to 



tr(BB^) 



(46) 



where r]k is the maximum allowed MSE at fcth user terminal. 
This problem can be written as the following optimization 
problem: 



mm 

B,G,C,r 

subject to 



Dfehfc - gfel 



< 



Vk, 



(47) 

1 < fc < M, 



< r. 



where r is a slack variable. With the NBE model of imperfect 
CSI, the robust transceiver design with MSE constraints can 
be written as 



mm 

b,g,c,r 

subject to 



(48) 



\\T>khk-gkf + (jl\\ck\\^ <ijk, 

Vhfee7^fc, i<k<M, 



In the above problem, the true channel, unknown to the 
transmitter, may lie anywhere in the uncertainty region. The 
transceiver should be so designed that the constraints are met 
for all members of the uncertainty set, TZk- This again, in the 
present form, is a semi-infinite optimization problem. In the 
following, we present a transformation that makes the problem 
in (Hill tractable. 

The optimization problem in ( l48T l is not jointly convex in 
b, g, and c. But, for fixed c, it is convex in b and g, and 
vice versa. So, in order to solve this problem, we propose 
an alternating optimization algor ithnll wherein each iteration 
solves two sub-problems. The first sub-problem given below, 
involves the optimization over {b, g} for fixed c. 



mm 

b,g,r 



(49) 



subject to 



IblP 



gfe|P+CT^||c;,^||2 <^^^ 

Vhfe e 7^fe, 1 < fc < M, 



< r. 



The second sub-problem involves optimization over {c} for 
fixed {b, g}, as given below 



mm 

c,si,--- ,sm 

subject to 



Sk 



(50) 



Vhfe eUk: l<k< M, 



where si, • • • ,sm are slack variables. The first sub-problem 
can be expressed as a semi-definite program (SDP), which is 

-For the case of single antenna users (i.e., MISO), a solution based on 
non-alternating approach is presented in [19]. 
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a convex optimization problem that can be solved efficiently 
[23]. Towards this end, we reformulate the problem in f4% as 
the following SDP; 



mm r 

b.g,r 



subject to 



Vk 

Dfc(hfc +efc) - gfc 

||b|| < r, 



Dfc(hfc +efc) -gfe 



(51) 



h 0, 



V|[efc|| < (5fc, 1 < < M, 



where r is a slack variable. In the reformulation given above, 
we have transformed the first constraint in (|49] l into an LMI 
using the Schur Complement Lemma [26]. We can show that 
the LMI in ( |5T] i is equivalent to 



A h P^XQ + Q^X^P, 



(52) 



where 



A = 



Dfchfc - gfc 

fnCfc 



Dfchfc - gfe 
cr„Ck 

I 



P = [O ] , X = Gfe, Q = - [1 0], and = 



Dfc 




Application of Lemma 1 to (|52] | and ( BTl i. as in Sec. IIV-AI leads 
to the following SDP formulation of the first sub-problem: 



mm r 

B,G,r,/3 



subject to 



rik - Pk 
Dfchfe - gfe 

0"nCfe 





Ibll <r. 



(Dfehfe - gfe 

OnCk 



I 



-5kTu 



(53) 



>- 0. 



Pkl J 

/3fe > 0, Vfc, 



Following a similar approach, it is easy to see that the 
second sub-problem can be formulated as the following convex 
optimization program: 



in b"+i and g"+i such that /fe(b"+\ g;!+\ c^) < r^k, 1 < 
k < M, where 



fk 



max Ek ■ 

hfcihfcGTJfc 



(55) 



Also, the transmit power P^+i = ||b"+i||2 < ||b"H2. Solving 
the second sub-problem in the n+lth iteration, we obtain c"+^ 
such that 



< A.(b"+\gr\<) 



(56) 



Since the transmit power Pr is lower-bounded and mono- 
tonically decreasing, we conclude that the sequence {P^} 
converges to a limit as the iteration proceeds. 

C. Robust MSE-Balancing Transceiver Design 

We next consider the problem of MSE-balancing under a 
constraint on the total BS transmit power in the presence 
of CSI imperfections. When the CSI is known perfectly, the 
problem of MSE-balancing can be written as 



mm 

B,G,C 



max efc 

k 



(57) 



subject to tr(BB^) < P„ 



This problem is related to the SINR-balancing problem due 
to the inverse relationship that exists between the MSB and 
SINR. The MSE-balancing problem in the context of MISO 
downlink with perfect CSI has been addressed in [30], [33]. 
Here, we consider the MSE-balancing problem in a multiuser 
MIMO downlink with THP in the presence of CSI errors. 
When the CSI is imperfect with NBE model, this problem 
can be written as the following convex optimization problem 
with infinite constraints: 



mm r 

b,g,c,r 



subject to 



Dfchfc - gfc 

l|b|| < 



(58) 



< r, Vhfc (^TLk, l<k< M, 



mm 

c, s,/^ 



Sfe 



subject to 



Sfe 



Dfehfe - 

0-„Cfe 





- gfe 



(Dfehfe - gfe 

0"„Cfe 



-SkT^ 





-SkFk 
Mfel . 



(54) 



>- 0, 



An iterative algorithm, as in Sec. IIV-BI which involves the 
solution of two sub-problems in each iteration can be adopted 
to solve the above problem. Transforming the first constraint 
into an LMI by Schur Complement Lemma, and then applying 
Lemma 1, we can see that the first sub-problem which involves 
optimization over b and g, for fixed c, is equivalent to the 
following convex optimization problem: 



mm r 



subject to 



The proposed robust MSE-constrained transceiver design al- 
gorithm alternates over both sub-problems. In the next sub- 
section, we show that this algorithm converges to a limit. 

1) Convergence: At the {n + l)th iteration, we compute 
b"+^ and g"+^ by solving the first sub-problem with fixed 
c". We assume that this sub-problem is feasible, otherwise the 
iteration terminates. The solution of this sub-problem results 



r~fJ,k 

Dfchfe - gfe 
cr,iCfe 




(Dfehfe - gfe 

CnCfe 



H 



I 



-(JfePfe 



1< /fc < M, 



(59) 



y 0, 



|b|| < VPn 



The second sub-problem which involves optimization over c, 
for fixed b and g can be reformulated as in (|50] |. By similar 
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Fig. 2. SMSE versus BS transmit power (Pt = ||B|||^) performance of 
tile proposed robust design in Sec. lIIII for ttie SE model. A't = 8, 6, M = 3, 

Nr^ = Nr2 = Nrs = 2, = L2 = L3 = 2, (T^ = ;^ = Q.l. 

Proposed robust design in Sec. 17/71 outperforms the non-robust design in [10]. 



Fig. 3. SMSE versus CSIT error variance (a-^) performance of tiie proposed 
robust design in Sec. |lIl]for the SE model. Aft = 6, 8, M = 3, = 
= Nr^ = 2, Li = L2 = Li = 2, Pmax = 15 dB, al = 1. Larger 
the value of cr ^, h igher is the performance improvement due to the proposed 
design in Sec. 17771 compared to the non-robust design in [10]. 



arguments as in the MSE-constrained problem, we can see that 
this iterative algorithm converges to a limit. 

V. Simulation Results 

In this section, we present the performance of the pro- 
posed robust THP transceiver algorithms, evaluated through 
simulations. We compare the performance of the proposed 
robust designs with those of the non-robust transceiver designs 
as well as robust linear transceiver designs reported in the 
recent literature. The channel is assumed to undergo flat 
Rayleigh fading, i.e., the elements of the channel matrices 
Hfc, 1 < fc < M, are assumed to be independent and 
identically distributed (i.i.d) complex Gaussian with zero mean 
and unit variance. The noise variables at each antenna of each 
user terminal are assumed to be zero-mean complex Gaussian. 
In all the simulations, all relevant matrices are initialized as 
unity matrices. The convergence threshold is set as 10^^. 

First, we consider the performance of the robust transceiver 
design presented in Sec. |III] for the stochastic CSIT error 
model. We consider a system with the BS transmitting L = 2 
data streams each to M — 3 users. In Fig. |2] we present 
the simulated SMSE performances of the proposed robust 
design and those of the non-robust design proposed in [10] for 
different numbers of transmit antennas at the BS and receive 
antennas at the user terminals. Specifically, we consider three 
configurations: i) Nt = 6, iV^ = 2, ii) Nt = 8, N,. = 2, 
and in) Nt = 8, TV^ = 3. We use cr| =0.1 in all the fliree 
configurations. From Fig. |2l it can be observed that, in all 
the three configurations, the proposed robust design clearly 
outperforms the non-robust design in [10]. Comparing the 
results for Nt = 6 and Nt = 8, we find that the difference 
between the non-robust design and the proposed robust design 
decreases when more transmit antennas are provided. A similar 
effect is observed for increase in number of receive antennas 
for fixed number of transmit antennas. It is also found that 
the difference between the performance of these algorithms 



increase as the SNR increases. This is observable in ( fTTI ). 
where the second term shows the effect of the CSIT error 
variance amplified by the transmit power In Fig. [3] we 
illustrate the SMSE performance as a function of different 
channel estimation error variances, (t|,, for similar system 
parameter settings as in Fig.|2] In Fig. [3] also, we observe that 
the proposed robust design performs better than the non-robust 
design in the presence of CSIT error; larger the estimation 
error variance, higher is the performance improvement due to 
robustification in the proposed algorithm. 

Next, we present the performance of the robust transceiver 
designs proposed in Sec. |IV] for the norm-bound model of 
CSIT error. Figure |4] shows the SMSE performance of the 
proposed design in Sec. IIV-AI as a function of the CSIT un- 
certainty size, 5, for the following system settings: Nt = 6, 4, 

M = 2, Nr^ = Nr^ = Nr^ = 2, Li = = ^3 = 2, 

Si — S2 — S, Pmax — 15 dB, and cr^ — 0.1. It is seen 
that the proposed design in Sec. IIV-AI is able to provide 
improved performance compared to the non-robust transceiver 
design in [10], and this improvement gets increasingly better 
for increasing values of the CSIT uncertainty size, 5. In 
Fig. |5] we illustrate the performance of the robust MSE- 
constrained design proposed in Sec. |IV-B] for the following set 
of system parameters: Nt = 4, 6, M — 2, Nr^ — Nr^ — 2, 
Li= L2 = 2, and 5i = 82 ^ 5 = 0.05, 0.1. We plot the total 
BS transmit power, Pt — ||B|||n, required to achieve a certain 
maximum allowed MSE at the user terminals, 7/1 — rj2 — rj. 
As expected, as the maximum allowed MSE is increased, the 
required total BS transmit power decreases. For comparison 
purposes, we have also shown the plots for the robust linear 
transceiver design presented in [20] for the same NBE model. 
It can be seen that the proposed THP transceiver design needs 
lesser total BS transmit power than the robust linear transceiver 
in [20] for a given maximum allowed MSE. The improvement 
in performance over robust linear transceiver is more when the 
maximum allowed MSE is small. 
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Fig. 4. S]y[SE versus CSIT uncertainty size (S) peii'ormance of the proposed 
robust design in Sec. ITV-Al for tlie NBE model. Nt = 6, 4, M = 2, = 

Nr^ = = 2, Li = L2 = L3 = 2, 5i = <52 = S, P^ax = 15 dB, 



0.1. Proposed robust design in Sec. \IV-A\ performs better than the 
design in [10]. 



Further, in Fig. |6l we present the total BS transmit power 
required to meet MSB constraints at user terminals for differ- 
ent values of CSIT uncertainty size Si ~ 82 = 5, for Nt ~ 4, 



M 



2, Nr 



Nr, = 2, 



Li = L2 = 2, and maximum 



allowed MSEs rji = ri2 = rj = 0.1,0.2,0.3. As can be seen 
from Fig. |6] the proposed robust THP transceiver design in 
Sec. lIV-Bl meets the desired MSB constraints with much less 
BS transmit power compared to the robust linear transceiver in 
[20]. We note that infeasibility of robust transceiver design for 
certain realizations of channels is an issue in robust designs 
with MSB constraints. In Fig. |7] we show the performance of 
the proposed THP transceiver design in Sec. IIV-BI in terms 
of the fraction of channel realizations for which the design 
is infeasible for Nt = 4, M = 2, Nr^ = Nr^ = 2, and 
different values of 6 and rj. It is seen that that the fraction of 
infeasible channel realizations in case of the proposed THP 
transceiver is much less compared to that in case of the linear 
transceiver in [20]. For example, for CSIT uncertainty size 
61=62^6 ^ 0.08 and user MSB rji ^ r]2 ^ V ^ 0.05, the 
linear transceiver design fails to produce a feasible solution 
in about 44% cases, whereas the proposed THP transceiver 
design fails only in about 24% cases. In Fig. [8] we show 
the convergence behavior of the proposed design. The number 
of iterations to converge depends on the MSB constraints. 
Stricter MSB constraints lead to larger number of iterations to 
converge. For example, when the MSB constraint is 77 = 0.3, 
the algorithm converges in around 6 iterations. Whereas, for 
1] = 0.1, it takes around 12 iterations to converge. 

Finally, in Fig. |9] we present the performance of the 
proposed robust MSB-balancing transceiver design in Sec. 
HV^ for Nt = 4, M = 2, Nr^ = Nr^ = 2, Li = L2 = 2, 
5i =62 = S = 0.02,0.1,0.15. The min-max MSB plots as a 
function of total BS transmit power constraint are shown. The 
corresponding plots for the robust linear transceiver in [20] are 
also shown. The results in Fig. |9] show that, for the same BS 
transmit power constraint, the proposed robust design in Sec. 



Fig. 5. Total BS transmit power (Py = ||B|||^) required as a function 
of maximum allowed MSE at the user terminals (rji = 172 = f?) in the 
proposed robust design in Sec. lIV-Bl for the NBE model. Nt = 4, 6, M = 2, 
= Nr^ = 2, Li = L2 = 2, CSIT uncertainty range 81=62 = 
S = 0.05, 0.1. Proposed robust THP transceiver design in Sec. \IV-B\ reauires 
lesser BS transmit power to meet the MSE constraints at the user terminals 
than the robust linear transceiver design in [20]. 



1<>- 











- Robust linear design in [20];ii=0.1 




8 




- e 


- Robust linear design In [20]; ii=0.2 








-0- 


Robust linear design In [20]:ti=0.3 




7 






- Robust THP design (sec.4.2) ri=0.1 


S 







- 


- Robsut THP design (sec.4.2) ii=0.2 


Q. 


6 




-0- 


Robust THP design (sec.4.2) r|=0.3 


1 










c 


5 








m 


4 








'b 





















H 


3 







— 

-e - 



- ^ 



-o- 



0.03 0.04 0.05 0.06 0.07 

CSIT uncertainty size, 8 



Fig. 6. Total BS transmit power (Pt = ||B|||^) required as a function 
of CSIT uncertainty size, 5, to meet MSE constraints in the proposed robust 
design in Sec. lR^ for the NBE model. Nt = 4, M = 2, N^ = Nr^ = 2, 
L\ = L2 = 2, maximum allowed JVISEs rji = ri2 = rj = 0.1,0.2,0.3, and 
<5i = ^2 = (5. Proposed robust THP transceiver design requires lesser BS 
transmit power than the robust linear transceiver design in [20]. 



IIV-CI achieves lower min-max MSB compared to the robust 
linear transceiver design in [20]. 

VI. Conclusions 

We proposed robust THP transceiver designs that jointly 
optimize the THP precoder and receiver filters in multiuser 
MIMO downlink in the presence of imperfect CSI at the 
transmitter. We considered these transceiver designs under 
SB and NBB models for CSIT errors. For the SB model, we 
proposed a minimum SMSB transceiver design. For the NBB 
model, we proposed three robust designs, namely, minimum 
SMSB design, MSB-constrained design, and MSB-balancing 
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Fig. 7. Fraction of infeasible channel realizations for different values of 
CSIT uncertainty size 5i = 52 = <5 in the proposed robust design in Sec. 
ITV-B] for the NBE model. Nt = 4, M = 2, = Nr^ = 2, Li = L2 = 2, 
and Tji = ri2 = r} = 0.05, 0.08. Proposed robust THP transceiver design 
in Sec. \IV-B\ has lesser infeasible channel realizations than the robust linear 
design in [20]. 



Fig. 9. Min-max MSE versus total BS transmit power limit, Pmax, in the 
proposed robust design in Sec. lTV-Cl for the NBE model. Nt = A, M = 2, 
Nr^ = = 2, Li = L2 = 2, (5i = &2 = S = 0.02, 0.1, 0.15. Proposed 
robust THP transceiver design in Sec. \IV-C\ performs better than the robust 
linear design in [20]. 
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Fig. 8. Convergence behavior of the proposed robust THP transceiver design 
in Sec. lTV-B] CSIT uncertainty size <5i = = <5 = 0.1. Nt = 4, M = 2, 
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2. Li = L2 = 2, and r)i = 7)2 = V = 0-1, 0.3. 



design. We presented iterative algorithms to solve these robust 
design problems. The iterative algorithms involved solution of 
sub-problems, which have either analytical solutions or can 
be formulated as convex optimization problems that can be 
solved efficiently. Through simulation results we illustrated the 
superior performance of the proposed robust designs compared 
to non-robust designs as well as robust linear transceiver 
designs that have been reported recently in the literature. 
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